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Abstract 

Basic aerodynamic coefficients are modeled as 
functions of angles of attack and sideslip with vehicle 
lateral symmetry and compressibility effects. Most of the 
aerodynamic parameters can be well-fitted using 
polymomial functions. In this paper a fast, reliable way of 
predicting aerodynamic coefficients is produced using a 
neural network. The training data for the neural network 
is derived from wind tunnel test and numencal 
simulations. The coefficients of lift, drag, pitching 
moment are expressed as a function of alpha (angle of 
attack) and Mach number. The results produced from 
preliminary neural network analysis are very good 

Introduction 

Wind tunnels use scaled models to characterize 
aerodynamic coefficients. The wind tunnel data, in original 
form, are unsuitable for use in piloted simulations because 
data obtained in different wind tunnels with different scale 
models of the same vehicle are not always consistent. 
Fitting a smooth function through the wind tunnel data 
results in smooth deri vativ es of the d^a. The smooth 
derivatives are important in performing stability analyses. 
Traditionally, the approach considered to descnbe the 
aerodynamics of the vehicle included developing 
wherever possible, a polynomial description of each 
aerodynamic function (B. Jackson and C. Cruz 1992). This 
ensured a smooth continuous function and removed some 
of the scatter in the wind tunnel data. Also, measurements 
of the same coefficient from two different wind tunnels are 
usually taken at dissimilar values of angle of attack and 
sideslip, and some means of reconciling the two dissimilar 
sets of raw data were needed. This curve fitting procedure 
is unnecessary for few coefficients. The curve fitting 
method used to generate the parameters for each 
polynomial description is an unweighted least squares 
algorithm. For the most part, tlie polynomial equations are 
generated using sparse data from wind tunnel experiments. 


Due to sparcity of data, mostly it will be defined as a line^ 
type of function. When more data are available, fhght 
control system designs will need to be revisited to allow 
for minor nonlinearities in control effects. 

Wind tunnel testing can be slow and costly due to high 
personnel overhead and intensive power utilization. 
Although manual curve fitting can be done, it is highly 
efficient to use a neural network (Magnus Norgaard, 
Jorgensen C and James Ross 1997, M.M.Rai and 
N.K.Madavan 2000 and Ching F.Lo, J.L.Zhao and 
DeLoach R 2000) to define complex relationship between 
variables. This paper is organized as follows: A short 
introduction to neural network and learmng followed by a 
section that will introduce the data set. Then results will be 
discussed to find an optimal solution to the various 
aerodynamic coefficients. The final section will conclu e 
optimizing the neural network and research directions. 

Neural Network 

A neural network is conceptually comprised of a collection 
of nodes and connectipn^(Rm5lh9Lf HmtM Wffiian^^^ 
1986; Wasserman 1989; Simpson 1990; Lau 1992; Masters 
1993; Jondarr 1996; Kartalopoulos 1996). The basic 
elements of a network are called neurons; tliey represent 
the sites that process information. The interconnecung 
links between the processing units, or neurons, are called 
synapses. Each synapse can be characterized by a weight, 
which is represented by a numerical value. All nemons in 
the adjacent layers are connected and the ° 

information is restricted to the forward direction. The 
network consists of three layers: input layer, hidden 
layer(s) and output layer. The hidden layer enables the 
network to learn relationships between input-output 
variables through suitable mappings. Among the many 
neural network models, the backpropagation algorithm is 
one of the better known and frequently used. 



Back propagation (Rumelhart, Hinton and Williams 1986) 
was created by generalizing the Widrow-Hoff learning rule 
(Widrovv and Hoff I960; Widrow and Stearn 1985; 
Widrow et al 1987) to multiple -layer networks and 
nonlinear differentiable transfer functions. Input - output 
pairs are used to train a network until it can approximate a 
function. Back propagation was the first practical method 
for training a multiple layer feed forward network. A 
neural network’s initial weights are simply random 
numbers, which change during training. Training consists 
of presenting actual experimental data to the neural 

network and using a mathematical algorithm the back 

propagation algorithm— to adjust the weights. Each pan of 
patterns goes through two stages of activation: a forward 
pass and a backward pass. The forward pass involves 
presenting a sample input to the network and letting the 
Ltivations flow until they reach the output layer. During 
the backward pass, the network’s actual output (from the 
forward pass) is compared with the target output and errors 
are computed for the output units. Adjustments of weights 
are based on the difference between the correct and 
computed outputs. Once each obseiwation’s computed and 
actual outputs are within the specified error tolerance, 
training stops and the neural network is ready for as-e; 
eiven a new input observation, it will estimate what the 
corresponding output values should be. After extensive 
training, the network eventually establishes the input- 
output relationships through the adjusted weights on t e 
network. 


Learning 


Learning in neural network is typically accomplished using 
examples. This is also called "traimng" of the neural 
network because the learning is achieved by adjusting the 
connection weights' in the neural network iteratively so 
that trained (or learned) neural network can perform certain 
tasks. Learning in neural network can roughly be divide 
into supervised, unsupervised, and reinforcement learning. 
Supervised learning is based on direct comparison between 
the acmal output of a neural network output and the desued 
correct output, also known as the target output. It is often 
formulated as the minimization of an error function such as 
the total mean square error between the actual output and 
the desired output summed over all available data, 
gradient descent-based optimization algorithm such as 
Lckpropagation can then be used to adjust connection 
weights in the neural network iteraUvely in order to 
minimize the error. Reinforcement learning is a special 
case of supervised learning where the exact desired output 
is unknown. It is based only on the information of whether 
or not the actual output is correct. Unsupervised leaiming 


is solely based on the correlations among input data. No 
information on "correct output" is avaikble for learmng. 
The essence of a learning algorithm is the learning mle, 
i.e., a weight-updating rule which determines how 
connection weights are changed. Examples of popul^ 
learning rules include the delta rule, the Hebbian rule the 
anti-Hebbian mle, and the competitive learmng rule. More 
detailed discussion of neural network and learning 
alsorithms are beyond the scope of this paper. The success 
of the neural network depends greatly on defimng the 
influencing parameters for the problem, definition o 
suitable neural network architecture and the learmng 
algorithm. 

Data Set for Aerodynamic Models 

Aerodynamic control systems can be divided into two 
categories viz,, control surfaces and aerodynaimcs 
controls. In this paper, aerodynamic controls and models 
are the focus. The variables involved in aerodynarmc 
controls are angle of attack (a), sideslip angle (p), elevon 
deflections (6e), aileron deflections (6a), rudder deflection 
(8R), speed brake deflection (5SB), landing gear effects 
and ground effects. The general equations of forces (lb) 
and moments (ft-lb) for key parameters are listed m the 
following tables 1 and 2 (B. Jackson and C. Cruz 1992). 


Forces (lb) 

Model 

Lift 

L=CL.q.S ^ 

Drag 

D=CD.q.S 

Side-force 

FY=CY.q.S 


Table 1 . Aerodynamic forces 


Moments (ft-lb) 

Model 

Pitching 

PM=Cm.q.S.c-i-(L.cosa-(-D.sinct).XMRC+ 

(L.sina-D.cosa).ZMRC 

Rolling 

RM-Cl.q.S.b+FY-ZMRc 

Yawing 

Cn. a. S.b-i-FY.XMRc 

■ 1 n X 1 . 


Table 2 Aerodynamic Moments 


The aerodynamic coefficients involved in the above 
equations are presented. 

Longitudinal aerodynamic coefficients 

Lift Coefficient CL; 

CL = CLBAs(a. M)-hACL.8H.APs(5F). 6 F-hACLspeedbra^(^ 
6SB)-HACLEG(5LG)+ACL,.(hfli)+ )+ACL,q(a, M).q.c/2U)+ 
ACL,„ (a, M).a’.c/2U 


Drag Coefficient CD: 


' Thresholds (biases) can be viewed as connection weights 
with fixed input 1. 





CD = CDbas(«. M)+ACD,8flaps(SF). 5F 

+ACDsPEEDBRAKE(a,8SB)+ACDLG(5LG)+ACD 


,e(h^)+ACD 


,q(a).q.c/2U 


Pitching Moment Coefficient Cm: 


Cm = CmB.\s(ot^ M)+ACm,gFLAPs(8F)' 

5F+ACmsPEEDBRAKE(^i 

8SB)+ACmLG(8LG)+ACmge(h/b))+ACm,q(a, M).q.c/2U)+ 

ACm,a(a. M).a’.c/2U 


Lateral aerodynamic coefficients 
Side force coefficient CY : 

(3Y-CYsb((^’ P+ACYsrudder (8R)- SR+ACYsaileron 

(5A)5A+ACYlga|j(5LG)P+ACYgea3(1'>^)P+'^^’^p(“)-P- 

b/2U+ACYXa).r.b/2U 

Rolling Moment Coefficient Cl: 

Cl=ClsB((“> - P+'^ClgRUDDER (8R)- SR+AClgAlLERON 

(5A)5A+AClLGAp(8LG)P+AClGEAp(h/b)P+AClp(a).p. 

b/2U+AClr(a).r.b/2U 

Yawing Moment Coefficent Cn. 

Cn=CnsB((f^^ P+ACn5RUDDER (SR)- 5R+ AC ns aileron 

(5A)5A+ACnLGAp(5LG)p+ACnGEAp(h/b)P+ACnp(a).p. 

b/2U+ACnr(a).r.b/2U 

Above equations depend basically on angle of attack and 
Mach number with little increments of other factors The 
above equation can be expressed as a function of angle of 
attack and Mach number and it resembles a simple 
polynomial expression. 


Depending on the geometry and mass properties of the 
vehicle, aerodynamics coefficients will vary. The genera 
parameters are tabulated in table 3. 


Parameters 

Ranges of values 

Angle of attack (degrees) 

-10<a<50 

Side angle (degrees) 

-20 < 3 < 20 

Mach number 

M<0.9 

Surface deflection (degrees) 

-15 < 5elevons (flaps) <15 
-20 < 6rudder < 20 
-20 < 5ailerons < 20 
0 < 5speedbrake < 80 

_ j font's 


Table 3 Range of values involved in aerodynamic coefficients 


If a regression analysis can fit a data set veiy well, then 
neural network can perform much better than regression 
techniques. Inputs considered for determimng base 


coefficients are angle of attack and Mach number. The 
output of the neural network is the coefficients o 
aerodynamic model. As a good training data set for a 
particular vehicle type, geometry and mass are selected 
from any wind tunnel test. Some times if the data set is not 
available from experiments for wind mnnels, a good 
training data set can be derived from numerical 
computations from Euler or Navier stokes or Vortex lattice 
method. This data set consists of a comprehensive mput 
and output tupple for an entire parameter space. Once 
training data set is defined, sparse data collected from 
experiments can be interpolated and extended for *e entiie 
range of data using a trained neural network. This wi 1 
avoid repeating the entire experiments in the wind ttinne . 
Once training data set is selected, one must determine the 
type of neural network architecture and transfer functions 
that will be used to interpolate the sparse data. The next 
section will discuss the selection procedure of the neural 
network architecture and transfer functions used in this 
work. 


N6ural Network Architecture 


ti this paper, interpolating for coefficient of lift is 
iscussed for sparse data set. The rest of the sanous 
erodynamic coefficients will be repeated with the same 
rchitecture of neural network with lespect to 
orresponding data set. The problem of defimng nemal 
etwork architectures (Freeman J. A and Skapma D.M 
992, Hagan M.T., Demuth H.B and Beale 1996) can be 
ivided into the following categories: (i) type of neiual 

etwork (whether three layer or four layer, etc.); (n) 
umber of hidden neurons; (iii) type of ti'ansfer functions 
Elliott D.L 1993); (iv) training algorithm; and (v) over 
nd under fitting of the results and validation of neinal 
letwork output. If the function consists of a finite number 
,f points, a three layer neural network is capable of 
earning the function. Since the availabihty of data is 
imited, the type of neural network considered for this 
ffoblem is a tliree layer neural network with input layer, 
tidden layer and output layer. The input layer mil have 
wo input neurons (alpha and Mach number) and output 
aver will represent one neuron (coefficient ot lit )• 
domain data has specific definite bounds, ratlier than 
laving no limits. The number of hidden neurons is defined 


For detennining an appropriate (hopefully optimal or netu- 
optimal) number of hidden units (Lawrence S, Lee Giles 
and Ah Chung Tsoi 1996), we construct a sequence of 
networks with increasing number of hidden neurons from 
to 20. More than 20 hidden neurons cause an over fitung 
of the results (Lawrence S, Lee Giles and Ah chung soi 
1997). Each neuron in the network is fully connected an 
uses all available input variables. First, a network with a 
small number of hidden units is framed using random 



initial weights. Iteratively, a larger network is constructed 
(up to the 20 hidden neurons) and the network results are 
compared with the expected results. Activation funcUons 
also play a key role in producing the best network results. 
The transfer function is a nonlinear function that when 
applied to the net input of a neuron, determines the output 
of the neuron. The majority of neural networks use a 
sigmoid function (S-shaped). A sigmoid funcuon is 
defined as a continuous real-valued function whose domain 
is the Reals, whose derivative is always positive, and 
whose range is bounded. In this aerodynamic problem, 
only a sigmoid function can alone produce an efficient fiu 
To «et a best fitting, it is suggested to use different kinds of 
tramfer functions for different layers of network 
However, functions such as “tanh” that produce both 
positive and negative values tend to yield faster training 
than functions that produce only positive values such as 
sigmoid, because of better numerical condiuomng 
Numerical condition affects tlie speed and accuracy of 
most numerical algorithms. Numerical condition is 
especially important in the study of neural network 
because ill-conditioning is a common cause of slow and 
inaccurate results from backprop-type algorithms. 
Activation functions for the hidden units are needed to 
introduce nonlinearity into the network. Without 
nonlinearity, hidden units would not make nets more 
powerful than just plain perceptrons (which do not have 
any hidden units, just input and output units). The reason 
is that a linear function of linear functions is again a linem 
function. However, it is the nonlinearity (i.e.. *e 
capability to represent nonlinear funcPons) that makes 
multilayer networks so powerful. Three types of acPvation 
functions are used in neural networks namely linear, 
sigmoid and hyperbolic tangent. 

The training epoch is restricted to 1000 cycles {present a 
data set, measure error, update weights}. The leammg rate 
and momentum is selected appropriately to get faster 
convergence of the network. The input and output values 
^e scaled to range [0.1, 0.9] to ensure that the output will 
lie in the output region of the nonlinear sigtiioid transfer 
function. Presentable variable values lie in between 0. 
and 0.9 (0.1 and 0.9 inclusive). The scaling is performed 
using the following equation 

A = r(V -V^) + 

X “ Amin 


r = 


V -V 

^ max nun 


Once scaled training data set is prepared, it is ready for 
neural network training. Levenberg-Marquwdt meAod 

(Martin T. Hagan and Mohammad B. Menhaj 1994)for 
solving the optimization is selected for back propagation 
training. It is selected due to its guaranteed convergence to 
a locarminima, and its numerical robustness. 

Experiments 

The training data set is divided into two sets viz,, dataset 
pairs which has Mach number less than 0.4, ^d those 
greater than 0.4^ The data set is presented to the neural 
Ltwork architecture for the training. Initially a training 
set which has 233 pairs, is presented to the neural network 
up ’to user-defined error of tolerance. The weights ^e 
stored and sparse data set of 9 pairs is provided for the 
same neural network architecture for further traimng. The 
initial training data set represents the exhaustive 
combination of data set in the particular parameter space. 
The initial training data set represents the general pattern 
of a particular aerodynamic coefficient. Based on the 
general pattern, the second training data set is interpolated. 
The initial data set is plotted in figure 1 and 2, and the data 
in fmure 1 can be represented by a linear type of funcUon 
whereas the data in figure 2 can be expressed as a 
combination of Unear and hyperbolic tangent or sigmoid 
function. Numerous trials have been conducted with 
different combinations of transfer functions, and we finally 
concluded that the linear transfer function be adopted or 
the input-to-hidden neurons and hyperbolic tangent or 
sigmoid function be used for the hidden-to-output layer. 
Figure 3 represents the sparse data set presented to the 
nemal network successively after the initial training a a 
set was presented. The figures 4 and 5 represent the neura 
network predicted data from the sparse data set. A few 
points are over fitted or under fitted m the results produced 
by the network. Over or under fitting is due to *e 
sparseness of data. Overall the results produced by the 
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Figure 1 Initial Training data 


for neural network (M < 0,4) 


V - Observed Variable 
A - Presentable Variable 


" Mach number < 0.4 then data expressed as a line^ 
function else a combination of linear and sigmoid function 




Conclusion 
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Figure 2 Initial Training data for neural network (M > 0.4) 
Alpha Vs CL 


0.8 
0.6 
d 0.4 
0.2 
0 


— Mac 0.3 
* ■ Mac 0.6 
--Mac 0.9 


8 

Alpha 


12 


Fi<nire 3. Sparse data presented to the neural network 
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Figure 4 Neural network interpolated data for sparse data 
(M < 0.4) 



Neural networks will become an important tool in future 
NASA Ames efforts to move directly from wind tunnel 
tests to virtual flight simulations. Many errors can be 
eliminated, and implementing a neural network can 
considerably reduce cost. Preliminary results have proved 
that neural network is an efficient tool to interpolate across 
sparse data. The prediction for the lower end and upper 
end of Mach number by the neural network is considerably 
deviated. The deviation is caused due to the non- 
availability of data in the sparse data. Imtially neural 
network has been trained by the original data which 
enables network to understand overall pattern. Successive 
training by the sparse data alters the weights of the neural 
network which causes this deviation. This deviation is 
well within 10 %, which is acceptable in aerodynamic 
modeling. 

Further reseai'ch is focused to overcome this deviation in 
predicting sparse data. It is also directed to optimize 
number of hidden neurons and will be integrated into web- 
enabled application. Hybrid system using evolutionary 
theory and neural network is planned to build an efficient 
model to predict aerodynamic variables. The neural 
network will be an integral tool of data mining suite in an 
existing collaborative system in NASA. 
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